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Abstract. Various sum rules accounting for the coupling between density and 
particle excitations and emphasizing in an explicit way the role of the Bose-Einstein 
condensation are discussed. Important consequences on the fluctuations of the par- 
ticle operator as well as on the structure of elementary excitations are reviewed. 
These include a recent generalization of the Hohenberg-Mermin- Wagner theorem 
holding at zero temperature. 
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1. INTRODUCTION 



The sum rule approach has been extensively employed in the literature in order 
to explore various dynamic features of quantum many body systems from a micro- 
scopic point of view (see for example [1] and references therein). An important 
merit of the method is its explict emphasis on the role of conservation laws and 
of the symmetries of the problem. Furthermore the explicit determination of sum 
rules is relatively easy and often requires only a limited knowledge of the system. 
Usually the sum rule approach is however employed without giving special emphasis 
on the possible occurrence of (spontaneously) broken symmetries. For example the 
most famous /-sum rule [2] holding for a large class of systems is not affected by 
the existence of an order parameter in the system. 

The purpose of this paper is to discuss a different class of sum rules which are 
directly affected by the presence of a broken symmetry. These sum rules can be 
used to predict significant properties of the system which are the consequence of 
the existence of an order parameter. In this work we will make explicit reference to 
Bose systems and to the consequences of Bose- Einstein condensation (BEC). Most 
of these results can be however generalized to discuss other systems exhibiting 
spontaneously broken symmetries. 

In sect. 2 we use a sum rule due to Bogoliubov in order to derive impor- 
tant constraints on the fluctuations of the particle operator as well as to obtain 
Goldstone-type bounds for the energy of elementary excitations. In sect. 3 we 
explore the consequences of BEC on the long range behavior of the half diagonal 
two-body density matrix and discuss the coupling between the density and particle 
pictures of elementary excitations in Bose superfluids. 
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2. THE BOGOLIUBOV SUM RULE. 



The Bogoliubov sum rule [3,4] 

/+00 
A-at,p('^)c^^ =< [flq - a-q,Pq] > = < ^0 + 4 >= 2y/ Nuq (1) 
-oo 

is the most remarkable among the sum rules depending explicitly on the order 

parameter tiq — \< gq >p. In eq.(l) AAt^si^) is the spectral function 

.4yit,B(w) = 4 V(e~^'^'" - e"^-^") < m \ \ n >< n \ B \ m > d{u - En - Em) 

m,n 

(2) 

relative to the operators A, B and we have made the choice = Cq — alq and S = 
Pq (Pg; and aLq are the usual density, particle annihilation and creation operators 
respectively and we have taken < ao >=< >= \/Nno). In deriving result (1) 
we have used the completeness relation and the Bose commutation relations for a 
and . Despite its simplicity, result (1) is at the basis of very important results 
characterizing the macroscopic behavior of Bose superfluids. Its richness is mainly 
due to the fact that it exploits in a microscopic way the effects of the "phase" 
operator, proportional to the difference Cq — alq. Some important applications of 
result (1) are discussed below. 

2.1. The Hohenberg-Mermin- Wagner (HMW) theorem. 

This theorem, of fundamental importance in condensed matter physics, states 
that no long range order can exist (no = in the Bose case) in a significant class 
of ID and 2D systems at finite temperature [5-7]. The theorem is based on the 
Bogoliubov inequality [3,4] 

< {A\A} >< [B\[H,B]] >> ksT \< [A\B] >|2 (3) 

3 



for the fluctuations of the operator A (we assume < A >=< B >= 0). Inequahty 
(3) can be derived starting from the inequahty [3,4] 



XAt,A<[B\[H,B]]>>\<[A\B]> 



(4) 



for the static response relative to the operator A 



XAt,A = 




AAi,B{'^)-duJ 



and the result < {A^ , ^} >> '2kBTxA^,A following from the fluctuation dissipation 



theorem. 



Choosing in the Bose case A''' = Oq and 5 = pq we obtain the inequality [5] 



n{q) =< a^Cq >> nom 



ksT 1 



(5) 



for the momentum distribution n{q) of the system. In order to derive result (5) we 
have used the /-sum rule [2] 



holding for Galilean invariant potentials. Ineq.(6) S{q,u!) = ^pt,p(<^)/(l— ea;p(— /?a;))| 
is the usual dynamic structure function. Inequality (5) points out the occurrence of 
an important infrared divergency in the momentum distribution which origi- 
nates from the thermal fluctuations of the phase of the condensate. This behavior 
permits to prove the absence of BEC in ID and 2D systems at finite temperature. 
In fact, due to such a divergency, the normalization condition for the momentum 
distribution X]qn(q') = N cannot be satisfied in ID and 2D, unless no = [5]. 

2.2. Extension of the Hohenberg-Mermin- Wagner theorem at T=0. 



u;S{q,u;)dw = - < [p^, [H,pq\] 




(6) 
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As clearly revealed by eqs.(3),(5) the HMW theorem does not apply at zero 
temperature. Actually there are important examples of 2D systems obeying the 
HMW theorem at T and exhibiting long range order in the ground state. The 
physical reason is that quantum fluctuations have in general a weaker effect with 
respect to thermal fluctuations in destroying the order parameter. 

An extension of the Hohenberg-Mermin- Wagner theorem holding at zero tem- 
perature has been recently formulated by Pitaevskii and Stringari [8]. To this 
purpose one uses the uncertainty principle inequality 

< {A\ A} >< {B\B} >>\< [A\B] >|2 (7) 

rather than the Bogoliubov inequality (3). Both inequalities provide rigorous con- 
straints on the fluctuations of the operator A. However, while the Bogoliubov 
inequality is sensitive to themal fluctuations and becomes less and less useful as 
T — > 0, the uncertainty principle is particularly powerful in the low temperature 
regime dominated by quantum fluctuations. By making the choice A^ = a^, B = pq 
one obtains the following non trivial inequality for the momentum distribution [8] : 

of a Bose system where S{q) — < PqPq > is the static structure function. The 
low q behavior of S{q) is fixed, at zero temperature, by the further bound 



S{q) < y J uS{q,u)du J —S{q,u)du! (9) 

that can be easily calculated at small q. In fact in this limit the inverse energy 
weighted sum rule (compressibility sum rule) J ^S{q,u!) = |x(?) approaches the 
compressibility parameter of the system (limq^oxio) = 1/mc^), while the energy 
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weighted sum rule (f-sum rule) is given by eq.(6). As a consequence S{q) < q/2mc 
at small q and the momentum distribution n(q) diverges at least as 

n{q) = no — . (10) 

It is worth pointing out that the 1/q law for the momentum distribution, already 
known in the literature [9,10], has been obtained here without any assumption on 
the nature of the elementary excitations of the system (phonons in a neutral Bose 
superfluid) and follows uniquely by the existence of BEC, the validity of the f-sum 
rule and the finitness of the compressibility. By imposing the proper normalization 
to the momentum distribution one can then rule out the existence of BEC in ID 
systems at zero temperature. Note that the ideal Bose Gas does not violate the 
theorem since in this case the compressibility sum rule x{q) diverges as 1/q^. 

Starting from the same inequality (7) it is possible to rule out the existence 
of long range order also in other ID systems at T=0 ( isotropic antiferromagnets, 
crystals) [8]. The theorem does not apply to ID ferromagnets since in this case the 
inverse energy weighted sum rule (magnetic susceptibility) diverges as 1/q^. 

Inequality (8) has been recently used [11] to rule out the existence of BEC in 
the bosonic representation of the electronic wave function 

*B(ri, ...,riv) = ea;p(- ^ai,j)^'F(ri, ...,rjv) (11) 

id 

proposed by Girvin and MacDonald [12] for the fractional quantum Hall effect. 
In eq.(ll) p is the fermionic wave function of electrons, v = ^^pi) where k is an 
integer, is the usual filling factor and aij is the angle between the vector connecting 
particles i and j and an arbitrary fixed axis. Using the Laughlin's expression [13] 
for the ground state wave function the authors of ref.[12] were able to show 
that there is not Bose-Einstein condensation in the bosonic wave function but 
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only algebraic long range order. The same result was obtained in ref.[14] starting 
directly from the Chern-Simons-Landau-Ginzburg theory. In ref.[ll] the absence 
of long range order was proven starting directly from the unceratinty principle 
inequality (8). The result follows from the fact that a charged liquid in an external 
magnetic field is characterized by a suppression of density fluctuations resulting in 
a quadratic law 



for the static structure function (uc is the cyclotronic frquency). This behaviour re- 
flects the incompressibility of the system and is consistent with the Kohn's theorem 
[15]. Inequality (8) implies a divergency for the momentum distribution of the 
bosonic wave function of this 2D problem. This is incompatible with the normal- 
ization of n{q), fixed by the total number of electrons, unless the Bose condensate 
relative to '^b identically vanishes. 

2.3. BEC and excitation energies. 

The sum rule technique can be used to to get useful constraints on the energy 
of elementary excitations. In particular use of the Bogoliubov sum rule (1) makes 
it possible to emphasize the role of the order parameter. 

Given a pair of operators A and B one can derive the following rigorous in- 
equality, holding at zero temperature, for the energy ujq of the lowest state excited 
by the operators A and B: 



Bounds of the form (13) were first considered by Wagner [4]. Result (13) can be 
obtained using the inequality (holding at zero temperature) 




(12) 




<[A\[H,A]] X [B^,[H,B]]> 
\< [A^,B] >|2 



(13) 




< [At, [if, A]] > 



(14) 



Xa^,a 
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where the r.h.s coincides with the ratio between the energy weighted and the in- 
verse energy weighted sum rules relative to the operator A. Use of the Bogoliubov 
inequality (4) then yields eq.(13). 

Result (13) has the important merit of providing a rigorous upper bound for coq 
in terms of quantities involving commutators. This is an advantage, for example, 
with respect to the so called Feynman bound: 

<[A^,[H,A]]> 
- <{^t,^}> ^1^) 

involving the anticommutator < {A^ , A} > in the denominator and hence requir- 
ing the direct knowledge of the fluctuations of the operator A. The occurrence 
of the anticommutator makes it difficult to exploit the low momentum regime of 
elementary excitations. Viceversa eq.(13) can be directly employed for this purpose. 

Result (13) is particularly interesting when the expectation value of the commu- 
tator is proportional to the order parameter of the problem. This happens 
in a Bose superfluid with the choice = Oq — alq and B — pq, already considered 
in this work. One then finds 

1 2 2 

ujl{q) < -|-[|- - + NV{Q) + y{p){n{^ + Q) + ^(P + Q))] (16) 
77-0 2,m Am ^ 

where we have carried out the commutators using the grand canonical hamiltonian 
B.' = H—jiN {jjL is the chemical potential) and taken a central potential with Fourier 
transform V{p). The quantitities n(p) and n{p) are defined by n{p) —< a|,ap > 
and n{p) = \ < f^pt^-p + OpO-p >• 

Result (16) is a rigorous inequality holding for any Bose system interacting 
with central potentials. It has the form of a Goldstone theorem. In fact, since 
its right hand side behaves as when — > 0, it proves the existence of gapless 
excitations, provided the order parameter uq is different from zero. 
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It is worth noting that in the hmit of a dilute Bose gas {fi = NV{0),n{p) = 
n(p) = NSpfi,no = 1), this upper bound coincides with the Bogoliubov dispersion 
law 

for any value of q. This result follows from the fact that the sum rules entering 
inequality (13) are exhausted by a single excitation, multi-particle states playing a 
negligeable role in a dilute gas. Of course in a strongly interacting system, such as 
liquid '^He, multiparticle excitations are much more important and the bound (16) 
turns out to be significantly higher than the lowest excitation energy u!o{q). 

It is useful to compare the Goldstone type inequality (16) with another, also 
rigorous, upper bound, still derivable from eq.(13), by choosing 

^ = q ■ jq - [if, pq] = ^ Y^iq' + 2q ■ k)al^^a^ (18) 

k 

and B = pq. In eq.(18) jq is the usual current operator. The resulting bound then 
coincides with the ratio between the cubic energy weighted and the energy weighted 
sum rules relative to the density operator pq [16,17] 

,.2(.^ <[[P-^,H],[H, [g,Pq]]]> 
< [p-q, [H: PqJJ > 

The explicit calculation of the triple commutator yields (we choose q along the 
^-axis) [16,17] 

< [[p_q,J/],[if, [ii,Pq]]]> 

= N^[{§^f + '^<Ek>+-^ [ ds{l - cosqz)g{s)VlV{s)] (20) 
m 2m m m J 

where < Ek > is the ground state kinetic energy and g{s) is the pair correlation 
function. The bound (19) then becomes: 

u^liq) < {^f + '^<Ek>+-^ f ds{l - cosqz)g{s)WlV{s) (21) 
Zm m m J 



and exhibits a quadratic behaviour in q as q ^ 0. It is worth noting that in the hmit 
of a dilute Bose gas (< Ek >= 0,^(s) = 1) both bounds (21) and (16) coincide 
with the Bogoliubov dispersion (17). 

The fact that it is possible to prove, via eq.(21), the existence of gapless excita- 
tions without assuming the existence of a broken symmetry is directly connected (see 
eq.(19)) with the conservation of the total current jq=o, holding in translationally 
invariant systems. This discussion also suggests that Goldstone type inequalities of 
the form (16) are particularly useful in sytems where the current is not conserved 
and where consequently only the existence of a spontaneously broken symmetry can 
ensure in a simple way the occurrence of gapless excitations. This is the case, for 
example, of spin excitations in magnetic systems (spin current is not conserved) or 
Bose systems with random external potentials. 

Let us discuss for example the effects of an external potential of the form (we 
take t/_k = U^) 

Vext = ^^P^ (22) 

on the sum rules discussed above. A first important result is that the Goldstone 
type upper bound (16) is not directly affected by the external field because of the 
exact commutation property 

[Oq - alq, [VexU Oq " O-q]] = . (23) 

Viceversa the cubic energy weighted sum rule for the density operator pq gets an 
extra contribution from the external force given by 

< [[p_q, H], [Kxt, [H, Pq]]] >= ^ < - E ^'^kPk > (24) 

m ^ 

At small q the new term provides the leading contribution to the triple commutator 
(20). Since the quantity < [p_q, pqj] > is not changed by the external force, 
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being still given by the /-sum rule (6), the upper bound (19) no longer vanishes 
with q. This different behaviour is due to the fact that the current jq=o is not 
conserved in the presence of the external field (22). This result also reveals that 
the relationship jq = q (a^ — a_q) between the current and the gradient of 
the phase operator, holding in a translationally invariant dilute Bose gas, is no 
longer valid in the presence of an external potential and gives rise to a normal (non 
superfluid) component of the density of the system at zero temperatrure. This is 
discussed, for example, in the Huang's contribution to this workshop. 

3. THE HALF DIAGONAL TWO-BODY DENSITY MATRIX. 

In this section we discuss another sum rule, also sensitive to the presence of 
Bose Einstein condensation, given by [18] 

/+00 2 
1 _ ^-pu: ^a+a^,p{^)d^ =< (Oq + a-q)Pq > (25) 

where A is the spectral function already introduced in sect.l with the choice = 
Oq + o^lq and B = pq. 

Differently from eq.(l), accounting for the commutation relation between the 
density and the phase operators, this sum rule cannot be expressed in terms of 
a commutator. Physically it accounts for the coupling between the density of the 
system and the modulus of the condensate, proportional to Oq + a^^q. An interesting 
feature of this sum rule is that it fixes the long range behavior of the half-diagonal 
two-body density matrix 

p(2)(ri,r2;r'i,r2) =< *^(ri)*t(r2)^r(r;)^r(r2) > (26) 

The occurrence of BEC is in fact not only relevant for the long range behavior of 
the one-body density matrix 

pW(r,r')=<*^(r)*(r') > (27) 
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fixed by the law 

Urn p(^)(r,r') =pno, (28) 

but also for the one of the two-body matrix (26). The LRO in the 2-body density 
matrix (26) is naturally defined by the equation [19] 

lim p(2)(ri,r2;r;,r2) = nop'(l + Fi(|ri-r2 I)) (29) 

and is characterized by the condensate fraction no as well as by the function -Fi(r). 
The properties of this function have been recently investigated in ref. [18,19]. The 
link between the sum rule (25) and the long range behavior of p*-^-* is fixed by the 
relation 

< Pq(aq + alq) >= ^/W^oil + 2Fi{q)) (30) 

where Fi{q) is the Fourier transform of -Fi(r). At low q the behavior of this function 
is fixed by the density dependence of the condensate according to law [18] 

q^o q 2nomc op 

The occurrence of LRO in the two-body density matrix implies the existence of 
a non trivial relation for the chemical potential that can be obtained starting from 
the following expression holding in systems exhibiting Bose-Einstein condensation: 

= E{N) - E{N -1)=<H> _ <^oHao> ^^^^ 

< alHao > 

Result (32) for the chemical potential follows from the property that the p = state 
in a Bose superfluid plays the role of a reservoir and that consenquently adding (or 
destroying) a particle in this state yields the equilibrium state relative to the N + 1 
{N — 1) system. Starting from eq.(32) one easily derives the equation: 

<al[H,ao]> 1 ^ t ^ T/^ ^ rr>o\ 
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Using results (1) and (30) one then finds the exact non trivial relationship [18] 

fi = pj drV{r){l + F^{r)) (34) 

relating the LRO function Fi{r) to the chemical potential of the system. Result 
(34) holds for any Bose system exhibiting condensation and interacting with central 
potentials. 

The LRO function Fi (q) plays a crucial role in the coupling between the density 

I F >= ^ pq I > (35) 



and particle 



I P >= ^^«-q I > (36) 

\/n{q) 

states that provide natural approximations to the elementary excitations of a Bose 
system. Result (35) coincides with the most famous Bijl-Feynman ansatz [20] . Both 
the density and particle pictures coincide with the exact eigenstates in the limit of a 
dilute Bose gas. For a strongly correlated liquid they provide only an approximate 
description. The coupling between the two states is different from zero because of 
occurrence of the Bose-Einstein condensation and is given by [1] : 



The coupling turns out to be complete when g — [in fact in this limit one has 
S{q) = q/2mc,n{q) = nQmc/2q and Fi{q) = — |] showing in an explict way that in 
the hydrodynamic limit the density and particle picture of elementary excitations of 
a Bose condensed system coincide. The coupling between the density and particle 
pictures is at the basis of fundamental properties exhibited by Bose superfluids (for 
a recent exhaustive discussion see ref.[21]). 
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It is finally interesting to compare the average excitations energies of the states 
(35) and (36). Both energies provide a rigorous upper bound to the energy of the 
lowest excited state of the system. The energy of the Feynman state (35) is given 
by the most famous law [20] 

<F\H\F> 

^^^^^ = <F|F> = ^''^ 
and is expressed in terms of static structure function. It is well known that in liquid 
^He eq.(38) provides the exact dispersion in the low q phonon regime, while it gives 
only a poor description at higher momenta. 

The energy of the particle state (36) takes instead the form [1] 

ep{q) = 1^-7^ ^ (39) 

where fx is the chemical potential and the quantity W{q) is given by 

W{q) = y"dridr;p(2)(ri,0;r;,0)e-^^('"'i-'"i)V(ri) (40) 

with p(^)(ri, 0; r'l, 0) defined in eq.(26). 

It is also useful to note that, due to relation (33-34) for the chemical potential, 
the particle energy (39) has no gap at q = and vanishes linearly with q [1] . 

It is possible to obtain an explicit expression for the average value of the particle 
energy (39) in momentum space: 

This average is sensitive to the values of ep{q) in the interval of momenta where 
the quantity q^n{q) has a significant weight. In superfiuid "^He this corresponds to 
the range q = 1 — 3A~^ including the maxon and roton region. Using the operator 
identity [22] - al[H, Op] =Ek + 2V the average (41) takes the form 

ep = fi - Ek > +2 < V >) (42) 
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where < Ek > and <V> are the kinetic energy and the potential energy relative 
to the ground state. At zero pressure, where fj, = j^{< Ek > + < V >), eq.(42) 
yields ep = — <V>/N = 21r^ 22K in superfluid ^He. It is instructive to 
compare the above value with the corresponding average of the Feynman energy: 

e-F = — • (43) 

Using microscopic estimates for S{q) and n{q) we find = 24 — 25K a value rather 
close to ep. 

The fact that the energies of the density and particle states (35-36) turn out 
to be comparable in the most relevant region q = 1 — 3A~^ reveals the importance 
of a careful microscopic investigation of the coupling (induced by BEC) between 
the two pictures in order to get a better understanding of the nature of elementary 
excitations in Bose superfluids. 
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